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Abstract 

We calculate differential cross sections for proton-nucleus elastic scattering by 
using a Glauber theory in the optical limit approximation and nucleon distributions 
that can be obtained in the framework of macroscopic nuclear models in a way 
dependent on the equation of state of uniform nuclear matter near the saturation 
density. We find that the peak angle calculated for unstable neutron-rich nuclei in 
the small momentum transfer regime increases as the parameter L characterizing the 
density dependence of the symmetry energy decreases. This is a feature associated 
with the L dependence of the predicted matter radii. 

Key words: Dense matter, Saturation, Unstable nuclei, Elastic scattering 
PACS: 21.65.+f, 21.10.Gv, 25.40.Cm 



Email address: keiiida@riken.go.jp (Kei lida). 
Preprint submitted to Elsevier Science 



4 February 2008 



The equation of state (EOS) of nuclear matter is essential to understanding 
of the saturation property of atomic nuclei [1], the structure of neutron stars 
[2], the mechanism of stellar collapse [3], and the dynamics of relativistic 
heavy- ion collisions [4]. Generally, the energy per nucleon of nuclear matter 
can be expanded around the saturation point of symmetric nuclear matter as 
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Here wo, no and Kq are the saturation energy, the saturation density and the 
incompressibility of symmetric nuclear matter, n is the nucleon density, and 
q; = 1 — 2x is the neutron excess. The parameters L and So characterize the 
density dependent symmetry energy coefficient S{n): So is the symmetry en- 
ergy coefficient at n = no, and L = 3no{dS/dn)n=no is the symmetry energy 
density derivative coefficient (hereafter referred to as the "density symmetry 
coefficient"). As shown in Ref. [1] by describing macroscopic nuclear proper- 
ties in a manner that is dependent on these EOS parameters, empirical data 
for masses and radii of stable nuclei can provide a strong constraint on the 
parameters Wq, no and 5*0, while leaving Ko and L uncertain. We remark that 
the isoscalar giant monopole resonance in nuclei (e.g., Ref. [5]) and caloric 
curves in nuclear collisions (e.g., Ref. [6]) can constrain Kq only in a way that 
is dependent on models for the effective nucleon-nucleon interaction. 

The incompressibility Ko and the density symmetry coefficient L control in 
which direction the saturation point moves on the density versus energy plane, 
as the neutron excess increases from zero. This feature can be found from the 
fact that up to second order in a, the saturation energy Wg and density are 
given by 

Ws = wo + Soo^ (2) 



and 
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The influence of Kq and L on neutron star structure can be significant [7] de- 
spite the fact that these parameters characterize the EOS near normal nuclear 
density and proton fraction, which are fairly small and large, respectively, as 
compared with the typical densities and proton fractions in the central region 
of the star. 

In our previous investigation [1] we pointed out the possibility that future 
systematic measurements of the matter radii of unstable neutron-rich nuclei 
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help derive the value of L. As an example, the matter radii of unstable neutron- 
rich nuclei such as Ni and even heavier element isotopes are expected to be 
deduced from future measurements of proton-nucleus elastic differential cross 
sections that will be performed by using a beam of such nuclei incident on a 
proton target and detecting scattered protons. This possibility is supported by 
our finding based on macroscopic nuclear models, which reasonably reproduce 
empirical data for masses and radii of stable nuclei and allow for uncertainties 
in the values of Kq and L, that the matter radii calculated for unstable nuclei 
depend appreciably on L, while being almost independent of Kq. However, ex- 
traction of the matter radii from experimental data for proton-nucleus elastic 
differential cross sections and for interaction cross sections is not straightfor- 
ward in the sense that it requires the approximate scattering theory [8,9]. It 
is thus instructive to examine how the cross sections themselves are related to 
the parameter L in a proper theoretical framework. 

In this paper we focus on proton-nucleus elastic scattering and, in the optical 
limit approximation of the Glauber multiple scattering model [10], obtain its 
angular distribution from the nucleon distributions in a nucleus calculated in 
Ref. [1] in a way dependent on L and Kq. For sufficiently high proton incident 
energies and small momentum tranfers to validate the Glauber model, we 
find that the angle of the scattering peak decreases with L more remarkably 
for larger neutron excess, while the peak height increases with Kq almost 
independently of neutron excess. We suggest the possibility that comparison 
of the calculations with experimental data for the peak angle may be useful 
for determination of L. 

We begin with macroscopic nuclear models used in this and a previous work 
[1]. We describe a spherical nucleus of proton number Z and mass number A 
within the framework of a simplified version of the extended Thomas-Fermi 
theory [11]. We first write the total energy of a nucleus as a function of the 
neutron and proton density distributions n„(r) and np{r) in the form 

E = Eb + Eg + Ec + {A- Z)mnC^ + ZmpC^ (4) 



where 

Eb^ j d^rn{r)w{nn{r),np{r)) (5) 



is the bulk energy with the energy per nucleon w{nn,np) of uniform nuclear 
matter. 
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is the gradient energy with adjustable constant Fq, 
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is the Coulomb energy, and m„ and mp are the neutron and proton rest masses. 
We express w as [11] 
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where 
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and 
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are the potential energy densities for symmetric nuclear matter and pure neu- 
tron matter. This expression for w is one of the simplest parametrization that 
reduces to Eq. (1) in the simultaneous limit of n — > no and a — > 0. We then 
set the nucleon distributions ni{r) [i = n,p) as 



nj(r) 



0, 
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and in the spirit of the Thomas- Fermi approximation minimize the total energy 
E with respect to -Rj, tj and nj"^ with the mass number A, the EOS parameters 
{no.wo, So, Kq, L) and the gradient coefficient Fq fixed. By calculating the 
charge number, mass excess and root-mean-square charge radius from the 
minimizing values of R^, ti and and fitting the results to the empirical 
values for stable nuclei (25 < ^4 < 245) on the smoothed beta stability line, 
we finally obtain no, Wq, 5*0 and Fo for various sets of L and Kq. We remark 
that as a result of this fitting, the parameters ai-62 become functions of Kq 
and L, while we fix the remaining parameter 63, which controls the EOS of 
matter for large neutron excess and high density, at 1.58632 fm^ throughout 
this fitting process. 



4 



The macroscopic nuclear models used here can describe gross nuclear prop- 
erties such as masses and root-mean-square radii in a manner that is dependent 
on the EOS parameters, L and Kq. Notably, as in Ref. [1], these models predict 
that the root-mean-square matter radii depend appreciably on L, while being 
almost independent of Kq. However, there are some limitations in the present 
macroscopic approach. First, this approach works well in the range of a < 0.3 
and A > 50, where a macroscopic view of the system is relevant. Second, the 
nuclear surface is not satisfactory in the present Thomas-Fermi-type theory, 
which tends to underestimate the surface diffuseness and does not allow for 
the tails of the nucleon distributions. As we shall see, such an underestimated 
diffuseness has consequence to calculations of proton-nucleus differential cross 
sections. Lastly, no shell and pairing effects are included. 

We proceed to calculate differential cross sections for proton elastic scatter- 
ing off nuclei using the Glauber theory in the optical limit approximation. The 
elastic differential cross section at given momentum transfer q and incident 
proton energy Tp can be written as (e.g., Ref. [12]) 



with the elastic scattering amplitude, 

ik 
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Here, b is the impact parameter, hk — ^J{Tp/c-\-mpC)'^ — {nipcy is the in- 
cident proton momentum, rj = Ze^/hv with the incident proton velocity 
V — hkc/{Tp/c + mpc) is the Sommerfeld parameter. 



^c(q) = ^ exp 
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is the amplitude of the Coulomb elastic scattering, which we approximate as 
a usual Rutherford scattering off a point charge, and 

«XN(b) = - y dv[np{v)Vpp{h - s) + n„(r)rp„(b - s)] (15) 



is the phase shift function with the projection s of the coordinate r on a plane 
perpendicular to the incident proton momentum and with the profile function 
FpTv of the free proton-nucleon {pN) scattering amplitude, for which we use a 
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simple parametrization, 



rpAr(b) 



1 — iapN 
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where cXpAr is the pN total cross section, a^jv = — Iiiirpjv(0)/Rerpjv(0), and 
(5pN is the slope parameter. 

In calculating the differential cross section, we use a numerical code based 
on the Monte Carlo integration for evaluations of the phase shift function 
(15), which can be applied to an arbitrary form of nucleon distributions. The 
essential input in such calculations is a set of the parameters, ctpAr, (5pN and 
(TpN, that characterize the profile function (16). Here the values of UpN-, PpN 
and (TpTv at given incident proton energy Tp are taken from Rcf. [13]. We finally 
obtain the differential cross section as a function of the CM. scattering angle, 
^c.m., from the nucleon distributions (11) in a nucleus determined as functions 
of the parameters L and Kq in the macroscopic nuclear models. 

We start with proton elastic scattering off stable nuclei. Figure 1 displays the 
elastic differential cross sections obtained for ^^Ni {Tp = 1047, 400 MeV), ^^"^Sn 
{Tp = 800 MeV), and ^ospb {Tp = 800 MeV). In each panel the calculations 
were performed at {L,Ko) =(5,230), (50,230), (80,230), (50,180), (50,360) 
in MeV. The dependence of the obtained differential cross sections on the 
EOS parameters is not appreciable. Note that the present Glauber model is 
valid for sufficiently large Tp to justify the optical limit approximation and for 
sufficiently small angles to allow us to ignore the nucleon-nucleon correlations 
in a nucleus. By comparison with the empirical data, we may conclude that the 
Glauber model reasonably describes proton-nucleus elastic scattering for high 
incident energies {Tp > 500 MeV) and for small scattering angles (^c.m. ^ 10 



Generally [17], the peak angles are related to the nuclear radius, while the 
peak heights are related to the diffuseness of the nuclear surface. Our nuclear 
model, on the other hand, predicts that the radius increases with L, while the 
diffuseness decreases with Kq. It is thus important to investigate the detailed 
peak structure in the small scattering angle regime and its relation with L 
and Kq within the present theoretical framework. 

In Fig. 2, the angles and heights of the first scattering peak, calculated for 
228 sets of L and Ko ranging < L < 175 MeV and 180 MeV < Kq < 360 
MeV, are plotted for proton elastic scattering off stable nuclei, ^^^Sn and ^^^Sn, 
at Tp — 800 MeV. ^ We first note that the peak angle decreases with L. This 

^ In this paper, we define the zeroth peak as that whose angle corresponds to 
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Fig. 1. Proton-nucleus elastic differential cross sections for ^^Ni (Tp = 1047, 
400 MeV), i^^Sn (Tp = 800 MeV) and ^o^Pb (Tp = 800 MeV). The solid, dot- 
ted, short-dashed, long-dashed and dot-dashed hnes are the results calculated 
for (L, Ko) = (50, 230), (50, 180), (50, 360), (5, 230), (80, 230) in MeV. For com- 
parison, we also plotted the results (thin lines) calculated from the Fermi-type 
nucleon distributions derived in Ref. [14]. The experimental data (dots) are 
taken from Refs. [14,15,16]. 



correlation is larger for larger neutron excess, since the matter radii increase 
with La^ [1]. Second, we find that the peak height increases with Kq in a 
way almost independent of neutron excess. We remark in passing that there is 
no appreciable correlation between the peak angle and as well as between 
the peak height and L. This absence of appreciable correlation between the 
peak angle and Kq is consistent with the fact that within the present Glauber 
model, the Fermi- type nucleon distributions derived by Ray et al. [14], which 
have larger surface diffuseness than those used in the present calculations while 
having similar root-mean-square neutron and proton radii, provide the angle 
of the first scattering peak similar to our results (see Fig. 1). We also note that 
Kq (L) induced uncertainties in the peak angle (height) are typically ±0.03 
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Fig. 2. The angles and heights of the scattering peak in the small angle regime, 
calculated as functions of L and Kq for p-^^^Sn and p-^^'^Sn elastic scattering 
at Tp — 800 MeV. The experimental angles and heights including errors (from 
Ref. [14]) are denoted by the horizontal hues (thick hues: central values, thin 
lines: upper and lower bounds). 



deg (±3 %), to which numerical errors due to the Monte Carlo integration are 
confined. 



The empirical data for the peak angle and height, also plotted in Fig. 2, 
contain uncertainties in the absolute scattering angle and normalization [14], 
which are here taken to be ±0.05 deg and ±10 %, respectively. The comparison 
between the calculated and experimental results for the peak angles suggests 
a possible extraction of L. Since the calculations ignore the pairing and shell 
effects and the tails of the nucleon distributions, however, the comparison of 
the absolute values of the peak angles is accompanied by systematic errors. 
These errors are expected to be reduced by considering a difference in the 
peak angle between the ^^^Sn and ^^^Sn cases. The calculated values of this 
difference depend on L only weakly and are within the uncertainties in the 
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Fig. 3. Same as Fig. 2 for p-^^Ni and p-^°Ni elastic scattering at Tp = 800 
McV. 

experimental value. We remark that extraction of Kq from the peak height is 
difficult, since the present calculations underestimate the surface diffuseness 
and the values of ctpjv, which control the peak-to- valley ratios [14], have yet 
to be determined uniquely. 

We turn to unstable nuclei, whose beams incident on proton targets may 
provide elastic scattering data in future experiments. In order to examine the 
possibility of extracting L from the peak angles for unstable nuclei such as ^°Ni 
(a = 0.3), we repeat the calculations for p-^^Ni and p-^°Ni elastic scattering 
at Tp = 800 MeV; the results for the angles and heights of the first scattering 
peak are exhibited in Fig. 3. We find that the L dependence of the difference 
in the peak angle between the ^^Ni and ^°Ni cases is appreciably large. The L 
induced change in the angle difference, AOc.m., amounts to order 0.1 deg. This 
is related to the L induced change in the matter radius difference, AR^^ ~ 0.1 
fm [1], by the diffraction pattern as /cA^^^.m. ~ t^AR^JR^^. In order to estimate 
L within ±20 MeV in the present theoretical framework, however, systematic 
measurements of the peak angles in the small momentum transfer regime for 
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various elements are desired for neutron excesses a < 0.3 with accuracy of 
±0.01 deg. 

In such an estimate of L, in which we focus on the angle of the first scat- 
tering peak, it would be useful to obtain a general relation between the peak 
angle and the parameters characterizing the adopted nucleon distributions 
within the Glauber model, as obtained by Amado ct al. [17] for intermediate 
momentum transfers and for a Fermi distribution. It is also important to ex- 
amine the predictability of the peak angle by the present scattering model. 
Neglect of the spin dependent part of the pN scattering amplitude and of 
nucleon correlations, as well as uncertainties in the values of ctpAr, may affect 
the prediction of the peak angle and hence the estimate of L, although they 
affect the prediction of the heights of diffraction minima and maxima more 
remarkably [18]. 

We finally consider proton-nucleus total reaction cross sections, which can 
be calculated from the same Glauber model as used for the elastic scattering 
calculations. The total reaction cross section can be written as 



where xn is given by Eq. (15). For neutron-rich stable nuclei such as ^osp]^ ^-^^^ 
Tp ^ 500 MeV, we find that the total reaction cross section tends to increase 
with either increasing L or decreasing Kq. This tendency stems from the fact 
that the reaction cross section becomes larger for larger diffuseness and radius 
of the nucleus. However, we have difficulty in deriving the EOS parameters 
from comparison with the empirical reaction cross sections even if data for 
unstable nuclei are accumulated. This is partly because errors in the empirical 
data [13] are of order or even larger than the change in the calculated reaction 
cross sections induced by the uncertainties in L and Kq and partly because 
the surface diffuseness and the tails of the nucleon distributions, which are 
underestimated and ignored in the present macroscopic nuclear model, affect 
the prediction of the reaction cross sections. 

In conclusion, by expressing the nucleon distributions in a way dependent 
on the EOS parameters and incorporating such distributions in the Glauber 
model, we calculated differential cross sections for proton-nucleus elastic scat- 
tering. At large neutron excess, we found out a large correlation between the 
peak angle in the small momentum transfer regime and the density symme- 
try coefficient L. This suggests a possible method to determine L from future 
systematic measurements of proton elastic scattering off unstable nuclei. 




(17) 
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